Cosmology with Hypervelocity Stars 
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In the standard cosmological model, the merger remnant of the Milky Way and Andromeda 
(Milkomeda) will be the only galaxy remaining within our event horizon once the Universe has aged 
by another factor of ten, ~ lO'^^ years after the Big Bang. After that time, the only extragalactic 
sources of light in the observable cosmic volume will be hypervelocity stars being ejected continuously 
from Milkomeda. Spectroscopic detection of the velocity-distance relation or the evolution in the 
Doppler shifts of these stars will allow a precise measurement of the vacuum mass density as well 
as the local matter distribution. Already in the near future, the next generation of large telescopes 
will allow photometric detection of individual stars out to the edge of the Local Group, and may 
target the ~ 10^*^ hypervelocity stars that originated in it as cosmological tracers. 

PACS numbers: 98.80.-k, 98.80.Es, 98.52.-b 



Introduction. In the standard LCDM cosmological 
model [H, the vacuum mass density which just started 
dominating the cosmic mass budget will remain constant 



= 0.7 X 10"^^ g cm 



(1) 



and so in the future, the scale factor of the Universe will 
grow exponentially with time a oc exp{_ffvt}, at a rate 
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(1.6 X 10^° yr)"\ (2) 



This accelerated expansion has important conse- 
quences for observers once the Universe ages by merely 
1-2 orders of magnitudes from now. Within ~ 10^^ yr 
after the Big Bang, all galaxies outside the Local Group 
will exit from our event horizon 0, 0] . At that time, the 
Local Group itself will be a single galaxy, Milkomeda, 
due to the imminent merger between the Milky Way and 
Andromeda in only ~ 5 x 10^ yr from now 4]. Milkomeda 
wiU be dominated by stars of low masses m^, ^ O.I-IM0, 
for which the lifetime t^ ?» 6.3 x 1O12(to^/O.1M0)-2-8 yr 
is comparable to the age of the Universe, t. Within 
^ 67 e-folds or ^ 10^^ yr, the wavelength of the cos- 
mic microwave background photons will be stretched by 
a factor of ^ 10^^ and exceed the scale of the horizon 
^hor '~ cll~^ — 4.9 X 10^ Mpc. At that time, not only ex- 
ternal galaxies but also the nearest extragalactic protons 
will be pushed out of our horizon and not be available for 
tracing the cosmic expansion. This realization has led to 
the naive expectation that an empirical reconstruction 
of the past cosmic history will become impossible 
Will future observers be unable to verify the validity of 
the standard cosmological model? 

Here we show that the continuous flow of hyperveloc- 
ity stars escaping Milkomeda will in fact allow future ob- 
servers to measure and the mass distribution within 
Milkomeda. This information will be sufficient to demon- 
strate the validity of the standard cosmological model at 
the above mentioned times. 



Cosmology with Hypervelocity Stars. Hypervelocity 
stars (HVSs) have velocities in excess of the escape ve- 
locity of their host galaxy. Such stars were discovered in 
the Milky Way halo over the past six years and are 
thought to be produced by stellar interactions with the 
nuclear black hole SgrA*. Although most of the known 
HVSs are bright and short-lived due to observational se- 
lection effects, one expects that lower mass stars are more 
commonly ejected from the Galactic center [lH. The 
ejection rate is low, once per ~ 10^^""^ yr, and is therefore 
expected to continue into the distant future long before 
the central reservoir of stars will be depleted in the Milky 
Way or Milkomeda. 

For a spherically-symmetric mass distribution, the or- 
bit of an HVS is purely radial, starting with some ejec- 
tion velocity > 10^ km s~^ at r ~ and subjected to two 
radial acceleration terms (l2| at r -C -Rhor, 



d r 



GM{< r) 



(3) 



with M{< r) being the galaxy mass interior to a phys- 
ical radius r. For a total Milkomeda mass of Mtot = 
10^^Mi3 Mq, the second ("cosmological constant") term 
will start to dominate at radii beyond 



-Rtran — 



/GMt, 



\ V 
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2.SM^l^ Mpc. 



(4) 



Figure 1 shows the relative contributions of the two accel- 
eration terms on the right- hand- side of equation ^ . The 
dominance of the second term at large distances moti- 
vated already the suggestion that the vacuum-dominated 
state of the future Universe could be probed with rocket 
flights [1^. Here we consider naturally produced HVSs 
in place of artificially made rockets to retrieve this infor- 
mation. 

With a typical escaping speed of ^ 10'^ km s^^, an 
HVS will traverse i?tran in a time ^ 2 x 10^ yr, much 
shorter than the lifetime of most of the stars ejected from 
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r (Mpc) 

FIG. 1: The net acceleration of an HVS in Milkomeda as 
a function of radius r (middle line). The Galactic deceler- 
ation (lower line) is calculated for a halo with a total mass 
Mtot = lO^^M©, a virial radius i?vir ~ 300 kpc, and a trun- 
cated density profile described by Eq. (1) of Ref. [l5|] (with 
i?200 = J?vir) that fits related numerical simulations. Half of 
the halo mass lies outside -Rvir- The acceleration due to the 
cosmological constant (upper line) dominates over the Galac- 
tic deceleration at large radii, r > Rtmn ~ 2.3 Mpc (see Eq. 



Milkomeda's center. By using advanced technologies, fu- 
ture observers could measure the acceleration of those 
stars from the change in the Doppler shift of their spec- 
tral lines. The distance of those stars can be found by 
inferring the intrinsic stellar luminosities from detailed 
spectroscopic observations. Measurements beyond i?tran 
will determine the value of p^, whereas closer HVSs will 
be used to find Milkomeda's mass distribution, M(< r). 

Since Milkomeda is surrounded by vacuum, it will have 
a sharp outer boundary 



15| . The boundary is interior to 



-Rtran bccause matter outside i?tran is accelerated out- 
wards and cannot be bound gravitationally. The aver- 
age matter density within the virial radius of Milkomeda 
Rvii ~ 300i?3oo kpc can then be measured. 



= 3 X 10^2^ 



M{< i?vir) 



5 X 1012 Mp 



^300 g cm 3. 



(5) 



The spherical collapse model implies that a galaxy em- 
bedded in a uniform matter-dominated background at 
early times should virialize with a mean density of ~ 100- 
200 times the background matter density, pm [l6[ . Future 
observers will then be able to estimate that the mean 
background density when Milkomeda had formed was 
(3 X 10-27 g cm-3/200) = 1.5 x lO'^s g cm'^. 
This value is coincidentally comparable to the vacuum 



mass density in equation ([T]). The age of the Universe 
when Milkomeda had formed can then be estimated to be 
^ H~^, since the total Hubble expansion rate H — a/a 
satisfies = (87rG'/3)(pni + Pv)- Because the matter 
density is diluted as pm oc a^^ , one will then be able to 
conclude that the Universe must have been dominated 
by matter at earlier times. 

The elapsed time since Milkomeda assembled most of 
its gas would be measurable from the age envelope of 
its stars and white dwarfs, t^. Such a measurement will 
provide a good estimate for the total time since the Big 
Bang, t « {H~^ + i*). The average matter density at t 
can then be verified to be unmeasurably small, Pm(t) ~ 
Pv exp{— 3(iJvt — 1)}. 

At t > 10^2 yr, the HVSs escaping from Milkomeda 
will be the only extragalactic sources of light filling up 
the observable volume of the Universe. For stars ejected 
with an initial speed from Milkomeda's center at r = 0, 
equation ([3]) can be integrated to give the velocity v{r) = 
(dr/dt) as a function of radius at distances r ^ i?hor. 
Given a steady ejection rate N, the radial profile of the 
number density of stars n(r) can be derived from the 
continuity equation. 



N 



47rr2ti(r) 



(6) 



At 
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^ ^tran the velocity of an HVS is reduced to 
(wq — w2g^)i/2^ where Vcsc is the escape speed from 
the production site of HVSs in the core of Milkomeda. 
Beyond i?tran the HVS velocity, v ^ [vtran + Hl{r'^ — 
i?2j.3^jj)]i/2, remains nearly constant out to a radius i?voi = 
(^'tranZ-H^v) = 16.3(i;tran/10'^ km s"^) Mpc, at which the 
cosmic acceleration starts to increase its value signifi- 
cantly. Therefore, at intermediate radii i?tran l^f^ Rvci 
equation ([6]) yields a power-law density profile, n oc r~^. 
This profile extends down to r = if ^ ^^csc- In 
this regime, the total number of HVSs increases lin- 
early with distance. However, at much larger distances 
-Rvci <^ i?hor the initial velocity is unimportant, and 
V sa H^r^ allowing one to measure pv from the slope of 
the Hubble diagram [v oc r) rather than from the accel- 
eration in equation ([3]). In this remote region n oc , 
and the total number of HVSs grows only logarithmi- 
cally with distance, providing a limited statistical ben- 
efit from improvements in the flux detection threshold 
of HVSs. Closer to the horizon, the cosmological red- 
shift (which is ignored here for simplicity) makes the 
total number counts even flatter with decreasing flux. 
For stellar masses > 0.5Mq, the HVS counts satu- 
rate at a sub-horizon distance for which the travel time 
equals the HVS lifetime. For example, solar-mass stars 
with a lifetime = lO^" yr can only shine while they 
travel out to a distance of ^ 30 Mpc for an ejection 
speed = 3 X 10"^ km s^^ and merely ^ 7 Mpc for 
= 10^ km s"!. 
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Since the observed energy flux from a star of a given 
luminosity is = L^,/(47rr^), the total number of 
HVSs in the sky above an observed flux F is given by 

N{> F)=4Tr nr^dr, (7) 

Jo 

where rniax(^") = (i*/47ri^)^/^. Figure 2 shows the re- 
sulting flux distribution, normalized by an HVS ejection 
rate of iV = 10~^ yr"^. For stars with > O.SM©, this 
distribution should be trimmed to flatten at the mini- 
mum flux corresponding to the distance for which the 
travel time equals the HVS lifetime. The travel time to 
a radius r can be read directly from Figure 2, since it 
equals (N/N) — lO^N yr in a steady state. Stars with 
= Q.5Mq and L^, = 0.07Lq (assumed to be forming 
steadily out of the gas reservoir in Milkomeda's nucleus) 
are the most luminous HVSs which could in principle be 
observed out to the horizon distance i?hor , since their life- 
time ti, ^ 7 X 10^" yr is comparable to the travel time 
from Milkomeda to ~ i?hor. A civilization living around 
an M-dwarf HVS with < 0.5Mq will have the privi- 
lege of taking a habitable one-way trip beyond the horizon 
of Milkomeda. They might find analogs of Milkomeda 
beyond its horizon, although they are most likely to find 
empty space since other galaxies are also separated away 
from them at an accelerated rate. 

Conclusions. In the distant future, runaway HVSs 
from Milkomeda will be the only extragalactic sources 
of light filling up the observable volume of the Universe. 
We have shown that a spectroscopic detection of the evo- 
lution in their Doppler shifts could be used by future ob- 
servers to validate the standard cosmological model even 
at a time when the wavelength of the relic radiation from 
the hot Big Bang exceeds the scale of the horizon. 

Indirect clues might also be available. The existence 
of an early radiation-dominated epoch could be inferred 
by measuring the abundance of light elements in metal- 
poor stars and interpreting it with a theory of Big Bang 
nucleosynthesis. The mass fraction of baryons within 
Milkomeda could be assumed to be representative of the 
mean cosmic value at early times. The nucleosynthesis 
theory can then be used to find the necessary radiation 
temperature oc a~^, such that the correct light ele- 
ment abundances would be produced. This would lead 
to an estimate of the time when matter and radiation 
had the same energy densities. Since density perturba- 
tions grew mainly after that time, it will be possible to 
estimate the amplitude of the initial density fiuctuation 
on the mass scale of Mtot that was required for mak- 
ing Milkomeda at a time {t — ti,) ^ H^^ after the Big 
Bang. Without a radiation-dominated epoch, this ampli- 
tude could have been arbitrarily low at arbitrarily early 
times. Future astronomers may already have cosmology 
texts available to them, but even if they do not, we have 
outlined a methodology by which they will be able to ar- 




logio(F/FMpc) = -21og,„(r/Mpc) 



FIG. 2: Number of HVSs out to a distance r (in Mpc) or with 
an observed flux > F (in unit of Fyipc = Li,/{A-R x Mpc^). 
Results are shown for an HVS ejection rate N = 10~^ yr~^ 
and the same mass distribution of Milkomeda that was used 
in Fig. 1. For long-lived stars with Li, — 0.071/0 (m* = 
0.5A/o and U ^ 7 x 10^° yr), the flux unit _Fmpc = 2.34 x 
10"^* erg cm^^ s^^ corresponds to a specific flux of ~ 0.8 njy 
or an AB magnitude of 31.7 at a wavelength of ~ 1/im and 
is comparable to the expected 1-a photometric sensitivity of 
the James Webb Space Telescope (JWST) [ij. The different 
lines correspond to different ejection speeds at r = 0, namely 
vo = 10^ km s~^ (upper line), 2 x 10^ km s"^ (middle hne), 
and 3 x 10^ km s~^ (lower line). 



rive at, and empirically verify, the standard cosmological 
model. 

The working assumption so far was that the vacuum 
energy density remains constant and does not decay to a 
lower energy state over hundreds of e- folding times. Mea- 
suring the dynamics of hypervelocity stars at late cosmic 
times provides the added benefit of testing whether the 
cosmological constant is truely constant. Indications to 
the contrary might reveal new physics, such as that con- 
sidered by evolving dark energy models [isj . 

One might wonder whether the signal shown in Fig- 
ure 1 can be used today to measure pv outside the Lo- 
cal Group. Conceptually, such a measurement would be 
contaminated by the large matter inhomogeneities that 
exist at distances smaller than tens of Mpc around the 
Virgo cluster. But even if one were interested in rnap- 
ping the mass distribution of those inhomogeneities [19j , 
the instrumental requirements are extremely challeng- 
ing. With state-of-the-art spectrographs using a laser 
frequency comb 12011 installed in the next generation of 
large telescopes [21| , one might hope to achieve at best a 
velocity precision of '-^ 1 cm s~^ , which corresponds to an 
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acceleration sensitivity of ~ 3 x 10~^ cm over a pe- 
riod of a decade. Unfortunately, this sensitivity threshold 
is still much larger than the expected acceleration ampli- 
tude of ^ 10~^° cm s~^ at the largest possible distance of 
~ 14 Mpc, which an HVS with a velocity of ^ 10'^ km s~^ 
could reach during the current age of the Universe. Al- 
though one could attempt to measure the velocity drift 
over a longer period of time, the actual acceleration sensi- 
tivity will be much worse than the best-case value stated 
above given the faintness of a single star at that distance 
(see Fig. 2). However, measuring the distance and ve- 
locity rather than acceleration of HVSs within the Local 
Group should already become feasible within the next 
decade. Future generations of experimentalists will have 

10^^ yr to improve upon this instrumental performance 
before it becomes crucial for validating the standard cos- 
mological model. 

The quantitative results presented in Figures 1 and 2 
may inspire related observational work in the near future. 
In particular, a ~ 2Mq star would be detectable photo- 
metrically by JWST (with a signal-to-noise ratio > 10) 
out to distances as large as a few Mpc. Thus, the next 
generation of large telescopes will allow photometric de- 
tection of individual stars beyond the edge of the Local 
Group and might target the ^ 10^^^ hypervelocity stars 
that originated in it as cosmological tracers. 
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